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SUMMARY 


Dynamical orbital techniques were employed to estimate 
the center-of-mass station coordinates of six C-Band radars 
located in the designated primary GEOS-C radar altimeter 
calibration area. This work was performed in support of 
the planned GEOS-C mission (December, 1974 launch). The 
sites included Bermuda, Grand Turk, Antigua, Wallops Island 
(Virginia), and Merritt Island (Florida). Two sites were 
estimated independently at Wallops Island yielding better 
than 40 cm relative height recovery, with better than 
10 cm and Im (relative) recovery for $ and A respectively. 

The tracking data used in this analysis were taken 
during 1969 when the radars tracked the GEOS- II transponder. 

The data used were exclusively that from the estimated 

sites and included 18 orbital arcs which were less than 

two orbital revolutions in length, having successive tracks 

over the area. In all, over 120 passes of data were 

used. Range biases were estimated. Error analysis and 

comparisons with other investigators indicate that better 

than 2m (1 a) relative recovery has been achieved at all sites. 

The techniques employed here, given their independence 
of global tracking support, can be effectively employed to 
improve various geodetic datums by providing very long and 
accurate baselines. C-Band data taken on GEOS-C should be 
employed to improve such geodetic datums as the European- 1950 
using similar techniques. 



INTRODUCTION 


The advent of artificial satellites permitted the science 
of geodesy to make measurements directly on a global scale. By 
studying satellite tracking data our knowledge of the Earth's 
geopotential has been increased enormously. In ever increasing 
numbers more and more precise tracking instruments are being 
deployed throughout the world in support of various planned 
satellite missions. Seemingly, the future continues to hold 
great promise for the science of satellite geodesy and geodyna- 
mics . 


The geodetic satellite missions planned for the remain- 
der of this decade will have the anticipated impact of greatly 
improving our knowledge of the Earth's size and geopotential 
while making valuable contributions to the fields of oceano- 
graphy and geophysics. The most immediate of these geodetic 
satellite missions is the Geodynamics Experimental Ocean Satel- 
lite: C (GEOS-C) with a planned December 1974 launch. The GEOS-C 
satellite will be extensively tracked by various metric sys- 
tems including C-Band and S-Band radars, lasers and various 
doppler instruments. An on-board radar altimeter will provide 
a means for measuring geoidal undulations represented by the 
sea surface's conformance to an equipotential surface. A 
global geoid of better than l®xl® resolution seems almost 
inevitable with this wealth of expected GEOS-C altimetry data. 

The SKYLAB altimeter data have already produced some important 
results . 

An important phase of the GEOS-C mission involves cali- 
bration of the altimeter system. This end will be partially 
achieved by studying the altimeter d'ta compared to a portion 
of the better known geoid. The western North Atlantic area 
has been selected for this purpose. A vast array of tracking 



instrumentation will be deployed at the tracking sites 
surrounding this calibration area for GEOS-C. Nevertheless 
in order to achieve satisfactory orbit determination over the 
calibration area, relative station positioning at the <2m 
level has become critical. Orbital error resulting from sta- 
tion positioning error even at this level will make the altime- 
ter calibration difficult. This report presents results which we 
believe satisfies the accuracy requirements for station posi- 
tioning around the GEOS-C calibration area. Figure 1 presents 
a map indicating the tracking sites of concern in this study 
and their location with respect to the calibration area. 



GODDARD 



FIGURE r 
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EXPERIMENTAL DESIGN 


During the months of January through February and 
later in October of 1969, the C-Band Network extensively 
tracked the C-Band beacon flown on the GEOS-II satellite. 

Table 1 presents the nominal orbit for GEOS-II, Incidental 
to this extensive radar tracking was some laser data taken 
at Goddard Space Flight Center (Station No. 7050), Mount Hopkins 
(7055) and Carnarvon (7054), Australia, A large amount of this 
radar data was taken by the sites f'f interest for the 
GEOS-C altimeter calibration. It was this data which were 
used to recover the center -of -mass positions of the radars 
located at Wallops Island (4840,4860), Bermuda (^.760), 

Grand Turk (4081), and Merritt Island (4082). Although 
Antigua (4061) has not been designated as a primary radar 
for the GEOS-C altimeter calibration, its geographical 
proximity to the calibration area caused us to select it 
for recovery also. Data from the Woomera (4946) C-Band 
site were also utilized for some related geopotential model 
analysis . 

Various techniques were considered for station position- 
ing. Dynamical orbital techniques were chosen primarily due 
to the suspected biases believed to be within the radar data. 

The GEODYN orbit determination system (T. Martin, 

1972) was used for the station recovery. GEODYN is a 
Bayesian least-squares, multiarc, multiple satellite orbit 
and geodetic parameter estimation system based upon Cowell 
type numerical integration techniques. Modeled parameters 
include luni-solar gravitational perturbations, solar 
radiation pressure, BIH polar motion and UTl data and 
several geopotential models . 
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TABLE 1. GEOS-II 


EPOCH 

APRIL 28, 1968 

Apogee Height 

1569 km 

Perigee Height 

1077 km 

Eccentricity 

O.CJ 

Inclination 

105.8 


Anomalistic Pariod 


112.1 min. 



Initial analysis centered upon two major concerns; 
selection of arc length and the best available geopotential 
model. Various arc lengths were simulated in an orbital error 
estimation scheme and it was found that an arc length of 
less than two revolutions had the desired effects of 
minimizing the force model error effects over the calibra- 
tion area when tracking was available on successive revolu- 
tions from the calibration area sites. In addition, no 
other data were required and even had detrimental impact when 
it was included into the orbital adjustment. This arc 
length is similar to the arc length Schutz, et al (1974) 
found to be optimal for their BE-C analysis. Very briefly, 
we found that this 1 1/3 revolution arc length and the 
selective inclusion of only calibration area data in the solu- 
tion had the desired effect of minimizing force model errors 
over the calibration area. By constraining the orbit in the 
same place on successive revolutions by including data in 
the solution only from the calibration area one can cause 
the force model error propagating with the frequency of 
the orbit to be minimized at the time of tracking. Figure 2 
presents a typical case of a GEOS- 1 1 orbit showing the geopoten- 
tial error propagation for an arc of the length described. 

Note that the dominant force model errors are of the period 
of the orbit and have been minimized over the calibration 
area for this arc length. 

The ORAN (Martin, 1970) error analysis system, which 
simulates a Bayesian least-squares adjustment, was used to 
create Figure 2. Twenty-five percent of two independently 
recovered geopotential models (APL 3.5 and SAO M-1; Martin 
and Roy, 1972) was propagated as representation of gravity 
model error. Figure 2 indicates that, as anticipated, the 
orbital error over the calibration area due to the geopoten- 
tial is only a very few meters. 



ORBITAL POSITION UNCERTAINTY (METERSI 


FIGURE 2. 


GEOS-II SATELLITE GEOPOTENTIAL ERROR PROPAGATION FOR TYrtCAL CASE 
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Thi*: arc length was also employed to evaluate various 
geopotential models. In a number of cases, ranging instru- 
ments located in Australia also tracked GEOS- 1 1 on the • 
revolutions as the calibration area tracking. The or’ . 
when passed through the Australian data unweighted in the 
solution, revealed the amplitude of geopotential error 
sources. Figure 2 shows that along track error due to geo- 
potential errors would be large over Australia. The along 
track error was measured by estimating the apparent timing 
errors in the residuals of the Australian data. 

Two state of the art gravity models were compared. 

They were : 

• The SAG Standard Earth II [SEII], Gaposchkin 
and Lambeck, 1970, and 

• The Goddard Earth Model 1 [GEMl] , Lerch, et al., 
1972. 

Table 2 presents a comparison of the recovered timing errors 
over Australia when these gravity models are used. GEMl was 
found to produce significantly smaller along track errors by 
almost a factor of two. Therefore, the GEMl geopotential 
model was initially adopted for our station recovery. As 
our work progressed, the GEM6 model (Lerci. , 197^) became 
available. Part of the Australian data analysis was repeated 
indicating comparable performance for this model with GEMl. 
However, since GEM6, unlike GEMl, used a large amount of 
metric data including GRARR, laser, doppler and C-Band obser- 
vations, it was this model which was employed for our final 
station recovery solution. 



TABLE 2. 

TIMING ERRORS OVER 
AUSTRALIA USING REFERENCE 
ORBITS DETERMINED WITH GEMl 
AND THE STANDARD EARTH II 
GRAVITY MODELS 

Epoch apparent 



Time GMT 

timing errors 

(msec) 

Station 

YYMMDD HHMM 

SE 11 

GEMl 

Woomera(4946) 

690201 2146 

-42.8 

-18.8 

Woomera(4946) 

690202 2206 

-6.9 

-2.1 

Woomera(4946) 

690203 0843 

14.6 

-0.9 

Woomera (4946) 

690204 2241 

6.9 

5.5 

Woomera(4946) 

69020S 0734 

3.4 

-0.7 

Woomera(4946) 

690205 0902 

10.7 

5.4 

Woomera(4946) 

690207 0957 

3.2 

2.9 

Carnarvon(7054 ) 

690207 0956 

7.3 

7.S 

Carnarvon(7054) 

690207 2342 

-62.2 

-34.4 

Carnarvon(7054) 

690208 1019 

7.2 

-2.7 

Woomera(4946) 

690211 2116 

-44.1 

-30.1 

Woomera(4946) 

690211 2306 

-14.3 

-12.0 

Woomera(4946) 

690212 0755 

-16.5 


Woomera(4946) 

690212 0945 

-11.8 


RSS 


93.86 

.59 
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DETAILS OF THE SOLUTION 


During most of the January/February period and all 
of the October period, the Eastern Test Range radars 
(Merritt, Grand '^ rk, and Antigua) tracked GEOS-II only 
during the daylight hours. Therefore, ot necessity, most 
of the arcs selected for our solution were from the daylight 
hours. This prevented us from including any laser data in 
our solution since only two daylight passes were available 
which were simultaneous with the selected best C-Band 
tracks and this was not enough data to recover the laser 
positions . 

We were fortunate however, for given the nodal rate 
of the GEOS-II satellite, the daylight passes taken in 
January/February were in the opposite direction from those 
in October (the later being South to North while Jan/Feb 
were North to South) . This provided us with favorable 
cancePation of the remaining geopotential error and good 
geometry for station recovery. Passes were selected on 
both sides of all the stations in both directions. 

Various data reduction experiments were performed to 
ensure that this data was corrected for known or suspected 
problems. The Wallops data required correction for timing 
errors of integer hundreds of milliseconds. The Bermuda 
observations from station 4740 were found to be inconsistent 
with the data from the other sites and since they were largely 
redundant with the data from 4760 they were deleted from the 
solution. The GODLAS laser data were found to have timing 
errors of about 10 msecs from February 5th through the I4th. 
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Eighteen arcs were selected for our solution. Table 
3 presents a list of the selected arcs. In all, over 120 
passes of radar data were used having a data sampling rate 
of 1 point/10 seconds. The final solution had over 6000 
range observations and had an RMS of fit of 1.4 m. Residual 
plots are presented in Appendix A. 
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TABLE 3. DATA ARCS SELECTED FOR C-BAND RADAR 
STATION ESTIMATION SOLUTION 


NUMBER OF RMS OF 

ARC NO . DATE AND HOUR OF EPOCH O BSERVATIONS FIT (M) 


1 

.'ANIJaRY 28, 1969 

9 HOURS 

308 

1.6 

2 

JANUARY 31. 1969 

8 HOURS 

270 

1.7 

3 

FEBRUARY 1, 1969 

8 HOURS 

212 

1.3 

4 

FEBRUARY 1, 1969 

22 HOURS 

349 

1.4 

5 

FEBRUARY 2. 1969 

9 HOURS 

297 

1.6 

6 

FEBRUARY 2, 19C9 

23 HOURS 

254 

1.5 

7 

FEBRUARY 4. 1969 

21 HOURS 

271 

1.9 

8 

FEBRUARY 8, 1969 

23 HOURS 

381 

1.9 

9 

FEBRUARY 12. 1969 

0 HOURS 

255 

1.4 

10 

FEBRUARY 12. 1969 

8 HOURS 

248 

1.9 

11 

SEPTEMBER 26. 1969 

15 HOURS 

455 

1.3 

12 

OCTOBER 8. 1969 

16 HOURS 

342 

1.1 

13 

r\r>--rr\r%Cr\ n 

UOI IOC 

3 in 

1 2 

14 

OCTOBER 10, 1969 

15 HOURS 

337 

1.4 

15 

OCTOBER 13, 1969 

16 HOURS 

366 

1.4 

16 

OCTOBER 16, 1969 

15 HOURS 

461 

1.2 

17 

OCTOBER 17, 1969 

16 HOURS 

525 

1.0 

18 

OCTOBER 24, 1969 

16 HOURS 

453 

1.3 


TOTAL 6095 1.4 
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RESULTS 


Table ^ presents the values for the C-Band station 
position recovery. The GEM6 gravity model was used for this 
final solution. The treatment of the C-Band biases are dis- 
cussed in the following paragraphs. 


The C-Band Biases 

The C-Band instruments have range biases which are 
significant when compared to their 1 to 1.5 meter noise 
level. A good general .nodel for a C-Band range observation 
would be : 


Rq » R + ARj + AR 2 + C 

where : 

Rq is the observed range 

R is the true range 

ARj is a bias term associated with very changeable 
errors so that 

ARj » Arj + * . . 

where are factors such as: 

Arj « errors due to thermal changes in the 
system 



TABLE 4 RECOVERED C-BAND STATION LOCATIONS 
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random calibration errors, etc. 


Ar, 


AR 2 is a bias term associated with long term ranging 
errors caused by unchanging systematic errors 
so that 


AR, : Ar 
where Ar 
Ar 

P 


P 

X 


♦ Ar . • . 

s 

are factors such as; 

a pulse width/bandwidth-mismatch problem 


Ar^ ■ an error to the surveyed ranging 

calibration target (or if the target 
is large, such as a water tank, an 
error to the tracking point) etc. 


£ is noise. If AR^ is small, the range measurements 
will be stable although they still may bo biased. 


Error analysis was performed and indicated that the 
station coordinates and all of the biases could be simul- 
taneously adjusted. An initial station recovery solution 
was made having all biases recovered independently for each 
individual pass of data. The recovered biases seemed very 
consistent from the Bermuda FPQ-0 (4760) and both Wallops 
radars. The systematic errors for these range passes exceeded 
the random errors at these sites making their observations 
very stable. We took advantage of this stability. The 
station recovery error analysis was again performed hut this 
time having the biases from Bermuda (4760) , and both Wallops 
stations (4840 and 4860) adjusted on an arc-by-arc basis. 



The results from the error analyses are presented in Table 5. 

The error sources modeled were: 

GMj^: 1 ppm error assumed 

Gravity: 254 of the difference between the 

APL 3.5 and SAO M-1 models. 

Timing: 0.1 msec error at all sites. 

Refraction: 54 error in tropospheric refraction 

(an error of 17 units of the refrac- 
tivity Ng) . 

The sites were adjusted with respect to Bermuda. Table 5 
indicates that superior results would be obtained if the 
Wallops and Bermuda biases could be adjusted on an arc by arc 
basis. This error analysis indicated that the 2m level of 
station recovery could be achieved. This station solution 
(which is quoted in Table 4) was then performed. 

Table 6 presents the actual differences in the recovered 
station positions when the bias treatment was modified as pre- 
viously discussed. The agreement in station recovery is consis- 
tent with the ORAN error analysis. We therefore were satisfied 
that our treatment of the C-Band biases did not adversely 
contaminate our station recovery. Our error analysis indicated 
that the 2m level of station positioning accuracy had been achieved 
at all sites. Our error analysis is probably conservative since 
the gravity model error propagated was scaled to the Standard 
Earth II and this model has been shown in Table 2 to be signifi- 
cantly less accurate (in the environment of this arc length) 
than either GEMl or GEM6. Gravity model error is still the 
dominant error source for this work. 

Table 7 presents the values of all the recovered biases 
from the adopted solution. 
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TABLE 5. ESTIMATl'O RELATIVE STATION RECOVERY UNCERTAINTY 

(METERS) 


BIASES FROM WALLOPS 




ALL BIASES 

ADJUSTED 

PASS-BYPASS 

AND BERMUDA ADJUSTED 
ARC-BY-ARC WITH ALL 
OTHERS PASS-BY PASS 

ANTIGUA 


1.2 

0.8 

(4061) 

A 

5.4 

1.9 


h 

2.2 

1.9 

GRAND TURK 

4> 

1.9 

o.e 

(4081) 

\ 

3.8 

2.1 


h 

2.0 

2.5 

MERRITT ISL. 


4.3 

1.3 

(4082) 

A 

1.6 

1.0 


h 

1.5 

2.4 

WALLOPS 


2.6 

0.7 

(4840) 

A 

2.9 

0.7 


h 

1.5 

1.6 

WALLOPS 

0 

2.7 

0.7 

(4860) 

A 

3.0 

0.6 


h 

1.0 

1.5 
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TABLP 6. DIFFERENCES IN RECOVERED STATION COORDINATES. 
SOLUTION WITH ALL BIASES ADJUSTED PASS-BY-PaSS MINUS 
SOLUTION WITH BERMUDA AND WALLOPS BIASES 
SOLVED FOR ON AN ARC BASIS 


STATION DIFFERENCE IN METERS 


NAME 

NUMBER 


__L. 

HEIGHT 

FTRANT 

4061 

0 15 

-0 79 

0.37 

ETRGRT 

4081 

0.22 

-0.69 

1.59 

ETRMRT 

4082 

-0.86 

-2.72 

-0.10 

NWALI8 

4840 

-0.61 

0.21 

1.28 

NWAU3 

4860 

0.06 

0.51 

-0.24 



TABLE 7 RANGE BIASES RECOVERED IN STATION bSTIMAllON SOLHTION 


STATION 


ARC NO. 

4061 

4081 

4^2 

4760 

4840 

4860 

1 

* 16 


8.7 

15.6* 


5 7* 


8.0 


-12.0 




2 

-36.6 



13 0* 

7.6* 



- 4.2 

5.7 





3 




16.0* 


11.0* 


-10.0 






4 

8.7 



18.0* 

4.8* 

6.2* 

5 

- 1.4 



20.0* 

5.7* 

7.S* 


- 0.0 






6 

- 6.8 



14.9* 

6.5* 

8.5* 

7 




18.7* 






5 7 



6.4 

8 

12.0 


8.4 

17.5* 


6.2* 




14.0 


12.0 


9 

0.8 


2.7 

12.4 

5.0* 





- 1.9 




10 

- 6.7 



8.8* 


1 

00 

• 


- 6.4 


-41.0 




11 

18.0 



18.0* 




- 0.0 

-14.0 





12 


4.0 

6.0 


- 0.3* 




- 0.7 

3.3 




13 


- 1.6 

4.1 



2.1* 




2.3 




14 

4.7 

11.9 

9.5 



- 7.7 


4.3 




- 1 6 


15 

11.2 

4.3 


- 8.5* 

- 0 6* 
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-11.1* 

1.4* 

1 7“ 


- 8.3 


- 6.2 




17 


- 3.9 


- ;.2* 


0 2* 



- 4.6 

0.0 




18 


11.8 

4.6 

1 

o 

• 

1.2“ 



3.3 - 0.8 

adjutttd on arc batii from two pattat of data. 





The stability of the Bermuda and Wallops radars are 
presented in Figures 4, and 5. In these figures the bias 
valies with their estimated uncertainties are presented for 
the successive arcs of data used in the solution except for 
arc 11 which is two weeks from any other data set. 

The bias recovery fjr the Wallops FPS-16 (Figure 3) 
deserves special mention. Farlv in 1969 a Frequency Shift 
Reflector (FSR) was installed at Wallops Island. Unfor- 
tunately the results from this point source target were 
unavailable during the January and February period. Therefore, 
this range data contained zero-set range errors due to cali- 
brating against a physically large reference range target. 

The magnitude of this error was determined by zero-setting 
the range with the range target and then calibrating this 
against the FSR. The apparent target size error found for 
the FPS-16 was 7,9m, This bias agrees extremely well with 
the values obtained from the station recovery solution 
(Figure 3). The October data set was corrected for this 
bias in the preprocessing. The bias recovered from the 
station solution had a mean of 0 meters during October. 

For all the radars it is not surprising to find that the 
mean bias values changed from Febiuary to October, since 
normal preventive maintenance could be expected to have 
this effect. 

Again, the stations were adjusted with respect to 
Bermuda. The error analysis indicated that any errors in 
the Bermuda position with respect to the center-of -mass 
would map virtually one to one into the recovered positions 
frem the other sites. Therefore, this did not effect our 
relative station po*^ itioning . In order to ensure center- 
of-mass recovery for the C-Band station coordinates, initial 
solutions were performed permitting Bermudans height to 
adjust . The longi tude and lat iiude were held unadjusted 
from center-of -mass values obtained 1 rom Marsh et al , 1973a. 



WALLOPS FPS-16 (4840) BIAS RECOVERY WITH ERROR 
ANALYSIS ESTIMATE OF UNCERTAINTY 
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The recovered height for Bermuda agreed to less than Im with 
the gravimetric geoid (Marsh, Vincent, 1973b). We wef" there- 
fore confident that we had a satisfactory Bermuda height as 
a reference point for our solution. Any error in the Bermuda 
longitude was not viewed as a problem since this merely es- 
tablished a longitude reference for tht calibration area. 
However, a center-of-mass latitude error at Bermuda would 
map into all stations nearly equally and therefore, would 
cause a systematic error in $ with respect to the center-of- 
mass. The accuracy quoted in the Marsh (op cit, 1973a) solution 
for Bermuda <p is 3m, and this uncertainty cannot be neglected. 
However, the relati/e station positi<^ning is not affected. 

Our error analysis indicated that it was desirable 
to hold at least one station height unadjusted in the solu- 
tion for this greatly lowered our sensitivity to variations 

in errors resulting from our adoption a priori of a value of 

3 2 

GMj;. A solution using a value of 39b601.2 Km /sec for GM^ 
and allowing Bermuda's height to adjust was compared with our 
final solution which held Bermuda's height fixed and used a 
value of 398600.8 Km^/sec^. After removing a difference of 
about 7 m in recovered height, the RMS agreement was 31, 24 
and 25 ,m in X, Y and Z, respectively. 


Results Compared To Ono-Pay Arc Solution 

We took an additional precaution to ensure center-of- 
msis recoverability. A solution using one- day arcs was also 
p,;rformed. Bermuda was again held unadjusted. The data which 
was selected was largely independent of the data used in the 
short arc recovery. The night-time GODLAS laser data was also 
included with the GODLAS position being adjusted. Sinc-j this 
data set was largely independent of the data used in the 
short arc recovery, we had a very small amount of data in 
this one-day arc solution from Grand Turk and Antigua. The 
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solution was therefore performed removing this Antigua and 
Grand Turk data since satisfactory station recovery from 
these sites was viewed as unlikely given their limited data 
sets. Table 8 presents the differences in X, Y and 2 between 
the short arc and one-day arc solutions. The recovered value 
for GODLAS is presented in Table 9 and is compared with the 
Marsh (op cit 1973a) recovery which used the same value of GM^ 
and was in the same longitude system. The agreement between 
the short arc and the one-day arc length recoveries is 
very good. We were also very satisfied with our recovered 
GODLAS position. 



TABLE 8. DIFFERENCES BETWEEN SHORT ARC AND ONE-DAY ARC 


SOLUTIONS IN METERS 




AX 

AY 

AZ 

STATION 

NUMBER 

im) 

itn) 


BERMUDA 

4760 

0. 

0. 

0 *h«ld untdjuttfd 

MERRITT 

4082 

-0.29 

-0.09 

-3.73 

WALLOPS 

4840 

0.05 

-0.30 

-1.11 

WALLOPS 

4860 

0.71 

-0.12 

-1.33 



TABLE 9. RECOVERED VALUES FOR GODLAS LASER STATION 7050 




LATITUDE 


LONGITUDE 


HEIGHT 


NUMBER 

DEG 

MIN 

SEC 

PEG 

MIN 

SEC 

(METCJ^jy 

ONE -DAY ARC 
RECOVERY 

7060 

39 

01 

14.183 

283 

10 

19.013 

4 15 

MARSH «t tl 
1973 

7060 

39 

01 

14.268 

283 

10 

18.955 

246 

DIFEERENCE 




-0/085 



0/058 

1 69 
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Comparison With Other Investigators 


In an investigation such as this one it is often very 
useful to comp«iie the station recovery results obtained uith 
those of some independent investigations. Extensive compari- 
sons were made to further assess the accuracy of the C*Band 
short arc recovery. The other investigators who were selected 
for comparison are presented below, along with a brief descrip- 
tion of their solutions. 

• Marsh, Douglas, Klosko, 1973: 

A global solution for stations using optical 
and laser measurements. 

• Icrch et al, 1974: 

GEM6 : a global solution using GRARR, laser, 

Doppler, C-Band and optical observations 
recovering station coordinates and the GEM6 
gravity model. 

• Martin and Walls, 1973: 

Narrow-band VLBI techniques were employed to 
recover values for the S-Band radars tracking 
the Lunar Excursion Module. 

• Walls, Boulware and Schanzle, 1973: 

S-Band station recovery from the tracking 
data taken on Mariner 9. 

• Marsh and Vincent, 1973: 

A gravimetric geoid combining satellite and 
surface measurements. 

• Ballew, 1974: 

A Gravimetric Geoid for GEOS-C Test Area based on 
satellite and surface measurements produced by the 
Defense Mapping Agency Aeronautical Center (DMAAC) . 



• Brooks 1972: 

A global solution using C-Band observrt io”s 
taken on GEOS-II for station coordinates. 

• Wolf Research and Development Corporation, 

1972 (Brooks): 

A solution in which the Grand Turk radar (4081) 
was recovered with respect to the Antigua radar 
(4061) held fixed at the Brooks (1972) value. 

The solution also recovered the geodetic position 
of a ship over the Irerto Rican Trench. 

• Air Force Eastern Test Ran-'e, 1974: 

Positions from the most recent "AFETR Geodetic 
Coordinates Manual." 

Figure 6 presents the geoid heights recovered by the 
various investigators. In order to account for the different 
values of GMg e.nployed, Figure 7 presents a plot of the geoid 
heights for the otner stations relative to the geoid height from 
Bermuda for each investigator. The gravimetric geoids agree very 
well with C-Band short arc values for geoid height except in the 
case of Bermuda relative to Wallops. In this latter case, all 
investigators who presented adjusted station coordinates 
are in substantial disagreement with tiie geoids as can be 
seen in Figure 7. The uncertainty quoted for the Marsh- 
Vircent gravimetric geoid is 2-4 m. The RMS agreement 
between the C-Band short arc solution and this gravimetric 
geoid is 2.8 m. The C-Band solution is even ir« closer 
agreement with the DMAAC geoid, especially at Antigua. 

Antigua is located in an area of large geoidal slope and 
DMAAC's use of IS'xlS' versus the Marsh-Vincent ' s l“xl* 
reduction may account for the difference seen. The RMS 
agreement between DMAAC and the C-Band solutions is 2.7m. 

The agreement between the Lerch solution and our solution 
is especially good. 
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HEIGHT RELATIVE TO BERMUDA (METERS) 


FIGURE 7. GEOID HEIGHT RELATIVE TO RERMUOA FROM VARIOUS 
INVESTIGATORS (BERMUDA'S GEOID HEIGHT MINUS OTHER SITE) 
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Table 10 presents the various recoveries for the 
Merritt Island station with respect to Bermuda. The com- 
pared quantities are relative spin axis distances and relative 
longitudes. In the cases of relative spin axis distances, the 
C-Band results show good agreement with the other investiga- 
tors. In the case of longitudes our agreement is best with 
the solutions with the smallest quoted uncertainties. Table 
ll compares the relative positioning of Grand Turk with respect 
to Antigua obtained from the Puerto Rican Trench Experiment 
and our C-Band short arc solution. 

Finally, the C-Band results have been compared to 
the North American Datum surveys. Since the Wallops radars 
were solved for independently*, a comparison can be made of 
their relative recovery with respect to their surveyed 
values. Table 12 presents these results. 

The North American Datum has recently been evaluated 
by using very precise traverses. A chord was measured up 
the East Coast of the United States. The C-Band recovered 
chord length is compared with this traverse in Table 13. 

Table 14 presents a comparison of the North American 
Datum chord from Wallops Island to Goddard Space Flight 
Center with the values recovered in the present analysis. 


*This comparison is especially interesting since 13 of the 
22 passes taken by 4860 were not simultaneous with data 
taken by the other Wallops radar (4840) . 
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TABLE 1C. RELATIVE SPIN AXIS DISTANCES AND RELATIVE LONGITUDES IMPLIED 
BY VARIOUS INVESTIGATORS BETWEEN BERMUDA AND MERRITT ISLAND 


INVESTIGATOR 

RELATIVE SPIN AXIS 
DISTANCE (BER-MER) 
(METERS) 

C-BAND RESULT - 
OTHER 

Marsh (opt/las) 

-219951.59 Hi 4m 

- 1.6 m 

Larch lopL ciop, las, CBND) 

-21 9955.08 6m 

1.9 m 

Martin (LEM) 

-219938.76 + 9m 

- 14.4 m 

Walls (Mariner 9) 

-219948.97 H_ 7m 

- 4.2 m 

Brooks 

-219957.23 + 7m 

4.1 m 

AFETR 

-219953.50 9m 

0.3 m 

C— Band Short Art 

-21 9953. 1 7 3m 




RELATIVE LONGITUDES 
IN DEGREES (BER-MER) 

C-BAND RESULT - 
OTHER 

Marsh 

16?0109778 4m 

-'!154 or —4.2 m 

1 erch 

16?01 09639 +. 6m 

104 or -2.8 m 

Martin 

16?01 10056 +. 7m 

~'!254 or -6.9 m 

Walls 

16?01 10211 i 7m 

-'!310 or -8.4 m 

Brooks 

16?011056 +. 7m 

-:'434 or -11.8 m 

AFETR 

ie?011043 *_ 9m 

-:'389 or -10 5 m 

C-Band Short Arc 

16?0109350 +. 3m 
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TABLE 11. RELATIVE LATITUDE, LONGITUDE AND HEIGHT FOR GRAND TURK 
AND ANTIGUA OBTAINED FROM THE PUERTO RICAN TRENCH EXPERIMENT 
AND THE C-BAND SHORT ARC RECOVERY 



RELATIVE 

RELATIVE 

RELATIVE 


LATITUDE 

LONGITUDE 

SPHEROID 


(GRT-ANT) 

(GRT-ANT) 

HEIGHT 

(GRT-ANT 

PUERTO RICAN 

TRENCH EXPERIMENT 

4.319000° 

-9.3395278° 

-7,96 m 

C-BAND SHORT ARC 

4.318947° 

-9.3395139° 

-9.29 m 

DIFFERENCE 

5.7 m 

-1.5 m 

1.34 m 
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TABLE 12. STATIONS 4840 AND 4860 RELATIVE RECOVERY 
COMPARED TO THE SURVEYED VALUES IN CENTIMETERS 


A A A Ah 

(cm I (cm) (cm) 


4840-4860 


-3 -73 -37 


TABLE 13. CHORD LENGTH COMPARISON BETWEEN 
MERRITT ISLAND AND WALLOPS 


CHORD PRECISE TRAVERSE C-BAND SHORT ARC 

4062 TO 4880 1149612.0 m 1149609.2 m 


TABLE 14. CHORD LENGTH COMPARISON BETWEEN 
WALLOPS ISLAND AND GODDARD SPACE FLIGHT CENTER 


CHORD NAD SURVEY 

7060 TO 4840 175728.61 m 


C-BAND SHORT ARC 
1/5730.21 m 
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2.8 m 
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-1.6 m 



CONCLUSIONS 


The station recovery presented in this paper is 
believed to satisfy the GEOS-C requirement of better than 
2m relative accuracy in each coordinate (la) for sites 
located in the GEOS-C altimeter calibration area. This 
conclusion is supported by our error analysis and com- 
parisons with other investigators. 

The techniques employed in this recovery experiment 
have wide-range implications for future GEOS-C experiments. 
The C-Band radars are all-weather instruments and unlike 
lasers, are not ad ersely affected by climatic conditions. 
Therefore* a geodetic parameter recovery experiment can be 
planned, analyzed, scheduled and the experimenter can be 
virtually certain that the necessary tracking data will be 
obtained by such radar systems. In addition, given that 
the techniques described within this report do not rely on 
global tracking support, a given experiment can be optimized 
v^^thout this difficult consideration. There will be a 
number of C-Band radars deployed on the European Datum for 
GEOS-C. The expected locations include Germany, France, 
Spain, Norway, and England. Radar data taken on GEOS-C 
should be similarly employed to recover accurate baselines 
for the reduction and improvement of the ED 1950. 
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APPENDIX A; RESIDUALS FROM C BAND STATION 
RECOVERY SOLUTION 



STATION KE'^ FOR RESIDUAL PLOTS 


SYMBOL 

STATION 

NAME 

LOCATION 

RADAR 

INSTRUMENT 

O 

ETRANT 

ANTIGUA, 

WEST INDIES 
ASSOCIATED STATES 

FPQ-6 

☆ 

ETRMRT 

MERRITT ISLAND, 
FLORIDA 

TPQ-18 

★ 

NBER05 

BERMUDA 

FPQ-6 

• 

ETRGRT 

GRAND TURK, 
BAHAMA ISLANDS 

TPQ-18 

■ 

NWALI8 

WALLOPS ISLAND, 
VIRGINIA 

FPS-16 

□ 

NWALI3 

WALLOPS ISLAND, 
VIRGINIA 

FPQ-6 
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